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In this article, we revisit the question of the valid- 
ity of Hartree-Fock and random-phase approximations. 
We show that there is a connection between the two and 
while the RPA as it is known in much of the physics lit- 
erature is of limited validity, there is a generalised sense 
in which the random phase approximation is of much 
wider applicability including to systems that do not pos- 
sess Fermi surfaces. The main conclusion is that the 
Hartree-Fock approximation is a mean-field idea applied 
to the density operator, and the random-phase approx- 
imation is a mean-field idea applied to the number op- 
erator. The generalised RPA is used to compute single- 
particle properties such as momentum distribution and 
spectral functions. It is found that we have to go beyond 
the generalised RPA and include fluctuations in the mo- 
mentum distribution in order to recover a nonzero imag- 
inary part of the one-particle self energy, which is also 
explicitly computed, all this in any number of spatial 
dimensions and no bosonization is needed. 



Both the Hartrcc-Fock Q| and the Random Phase 
approximations (RPA) 0] are widely used in the 
physics literature to study Fermi systems. The Bo- 
goliubov theory || is the analog of the random-phase 
approximation for Bose systems. This latter fact has 
been demonstrated in detail in our article Q]. The 
analysis presented here could be repeated for Bose 
systems as well. Here we try and address the ques- 
tion of validity of the Hartree-Fock approximation 
and the RPA. We show that the mean-field approx- 
imation carried out on the density operator is the 
Hartree-Fock approximation. The RPA manifests 
itself as mean-field theory applied not to the density 
operator but to the number operator. The density 
operator measures how the electrons are distributed 
in real space where as the number operator mea- 
sures how the electrons are distributed in momen- 
tum space. Just as the Hartree-Fock approximation 
is valid when fluctuations in the density of electrons 
at each point in real space is small, the RPA is valid 
when fluctuations in the momentum distribution of 
electrons is small compared with the average mo- 
mentum distribution which measures the probability 
of an electron to possess a given momentum. 

Let us try and write down some representative ex- 



amples when RPA and Hartree-Fock approximations 
are used. Let us take for example, the jellium model 
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If we apply the mean-field idea on the density, that 
is, replace p q by (p q ) then we get a hamiltonian that 
does not involve any coulomb interaction at all(apart 
from an additive constant). 
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Therefore this approximation is bad in the extreme. 
However the random-phase approximation is still 
valid for this system. It is the mean-field idea ap- 
plied not to the density but to the number operator. 
For this we have to rewrite the full hamiltonian given 
below, 

H = J2 e k4 c k + E W E C U q /2 C l'- q /2 C k'+q/2 C k-q/2 
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This has to be replaced by, 
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where 
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where n k = clc k - Let us write, 
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We plug the above decomposition into Hq and find 
that if we neglect terms quadratic in the fluctua- 
tions, we get what we are after, namely the gener- 
alised RPA called for simplicity as just RPA. 
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The average occupation is 
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The chemical potential fi has to be fixed by making 
sure that, 
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At zero temperature fi = ep, this quantity is equal 
to the usual Fermi energy when u q = 0. 



(n k ) = 0(e F - e k ) = 6>(e F - e k + 5] y (™k- q )) 
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and # is the Heaviside step function. We can 
now demonstrate that the generalised RPA dielec- 
tric function(the Lindhard dielectric function || be- 
ing a weak coupling limit of this) may be recov- 
ered using the following procedure. If one considers 
an extremely weak external perturbation applied to 
the system and follows the discussion in Mahan § 
one arrives at the following formula for the dielectric 
function, 



e fiPA (q,w) = 1 + y 5Z 
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The only point to bear in mind is that we have 
to use the full interacting momentum distribution 
rather than just the noninteracting value. The above 
formula differs from the traditional RPA in two re- 
spects. First, we have the full interacting momen- 
tum distribution in the numerator. This has to be 
determined self-consistently by an equation such as 
Eq.( |) or Eq.( While this feature may be found 
in our earlier work p|, the second feature is new. 
The denominator contains the renormalised disper- 
sion e k rather than the parabolic e k . Thus we can 
see that there is a whole new set of approximations 
that go beyond the RPA. While none of these revela- 
tions may come as a surprise to the reader, it should 
serve as a reminder that even our most cherished ap- 
proximations may not be controlled in any sense of 
the term. It is more likely that they were the first 
to appear in the literature and probably the easiest 



ones to use thus explaining their popularity. The 
dielectric function written down above has the at- 
tractive feature of reducing to the familiar Lindhard 
dielectric function [0 for extremely weak coupling 
and at the same time giving us something very dif- 
ferent for stronger coupling. The only drawback of 
the above approach is that if we compute the one- 
particle Green functions we find that the imaginary 
part vanishes identically. This is unfortunate and we 
have to do better in order capture lifetime effects. 
The RPA hamiltonian neglects fluctuations in the 
momentum distribution of the electrons. In order to 
recover a finite lifetime of single particle excitations, 
we find that it is important to study the generalised 
Hq rather than the more simple HppA ■ The fluctua- 
tions in the momentum distribution may be related 
to the mean by the following observation. Define, 
N(k, k ) = (n k n k <) — ("k)(«- k ') • The fluctuation 
in the number operator is N(k, k) = (n£) — (n k ) 2 . 
Since n 2 . = n k for fermions, we have 



iV(k,k) = (n k }(l-(n k » 
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Therefore, we may conclude that any nonideal mo- 
mentum distribution fluctuates (there are however, 
pathological exceptions, see footnote Q). In fact, a 
very nonideal momentum distribution such as one 
for which (n k ) = 0.5 for most momenta has the 
largest fluctuation. When dealing with nonideal sys- 
tems, we are obliged to consider fluctuations in the 
momentum distribution. In order to study lifetime 
effects therefore, we have to include the full Hq. This 
may be written more transparently as (apart from 
additive constants) Hq — HppA + Hfi, 
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The full Fermi propagator may be evaluated by 
treating the fluctuation part as a perturbation and 
using the functional methods of Schwinger illus- 
trated brilliantly by Kadanoff and Baym |7j] . Before 
we plunge into the details it is important to keep in 
mind that the mean also changes when we consider 
fluctuations. That is, 
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here (n k ) rpa is given by Eq. ( ||) or Eq. ( |ll|) the rest 
is nonzero only when the fluctuation in the momen- 
tum distribution is large (which is, unfortunately, 
almost always the case when (n k ) ppa is nonideal) 
Therefore, now, <5n k refers to fluctuation around the 
full average rather than the RPA average. The fi- 
nal answers are given below, and it is hoped that 
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the reader can rederive them using the references 
quoted in the bibliography (mainly || , Qj). We shall 
adhere to the notation of Kadanoff and Baym . In 
their notation the final answers for the single-particle 
Green functions are as follows(we assume in the fol- 
lowing that F (p) ^ 0, however one may investigate 
the limit F(p) -»■ 0, here z n = (2n + l)n//3), 
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From Eq. ( [l7|) we may obtain the real and imaginary 
parts of the retarded self-energy, and from there the 
spectral function and the collision rates . 



r(p,w) = tJ-k(p,u) 
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if k(p, uj) — (lu — e p + n) 2 — 4F(p) < and both are 
zero otherwise (when F(p) ^ 0). It can be seen that 
the spectral function is peaked around e p — /i with a 
width of the order of 2y/F(p), and the collision rate 
is vanishingly small for those values of p for which 
F(p) is close to zero. It may be shown that the mo- 
mentum distribution including possible fluctuations 
is given by, 
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Again, it may be seen quite easily that Eq.( [H]) is 
identical to Eq.( ||) when fluctuations in the mo- 
mentum distribution are ignored. The only sticking 
point now is the computation of the fluctuation in 
the momentum distribution. This may be done in 
a similar manner by employing the functional meth- 
ods of Kadanoff and Baym M. However, we shall 
instead adopt a simpler approach based upon some 
earlier work ||. There we showed that for k ^ k 
that the number fluctuation has this rather simple 
looking form, 



N(k,k) = -(n k )n p (k) 
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where n^(k) = l/(exp(/3(ek — /io)) + 1) and /io is 
the chemical potential of noninteracting fermions (at 
finite temperature). This leads to a rather compli- 
cated set of self-consistent equations for the momen- 
tum distribution and its fluctuation. These formu- 
las have very appealing and illuminating forms and 
they also support our earlier conclusions [Q namely 
that Fermi liquid theory can break down in all three 
spatial dimensions for sufficiently strong values of 
the coupling and be intact for sufficiently weak cou- 
plings. This analysis also provides another (apart 
from our earlier preprint [|j) illustration of the im- 
portance of the concept of fluctuations in the mo- 
mentum distribution. Lastly, it is worth remark- 
ing that we have not really answered the question 
posed in the title namely, is the random phase ap- 
proximation(RPA) a controlled approximation ? It 
is our belief that the answer is no, for RPA cannot be 
thought of as the leading behaviour of an expansion 
in powers of an obvious small parameter. The fact 
that the traditional RPA may be obtained as a weak 
coupling limit of Eq. ( [l2|) does not clarify the situ- 
ation, for it may well be that even the generalised 
RPA is not a controlled approximation. Indeed, we 
know for a fact that even the generalised RPA is not 
sufficient in giving us a finite lifetime. Thus the title 
is somewhat rhetorical and it is hoped that critics of 
our earlier work will take note. This work was sup- 
ported in part by ONR N00014-90-J-1267 and the 
Unversity of Illinois, Materials Research Laboratory 
under grant NSF/DMR-89-20539 and in part by the 
Dept. of Physics at University of Illinois at Urbana- 
Champaign. The authors may be contacted at the 
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1 See for example, " Quantum Mechanics" A. A. 
Sokolov, I.M. Ternov and V.Ch.Zhukovskii, Mir Pub- 
lishers, Moscow, 1984 

2 D. Bohm and D. Pines, Phys. Rev. 92, 609 (1953) 

3 Bogoliubov N. N., J. Phys. Moscow 11, 23 (1947), 
Bogolhibov N. N. and D. N. Zubarev, Sov. Phys. 
JETP 1, 83 (1955) 

4 G.S. Setlur and Y.C. Chang, Phys. Rev. B, 57, no. 
24, 15 144(1998) 

5 G.D.Mahan "Many-Particle Physics", II ed.,Plenum 
Press, New York, 1990 



3 



J. Lindhard, K. Dan. Vidensk. Selsk. Mat. Fys. 
Medd. 28, (8) (1954) 

L. Kadanoff and G. Baym " Quantum Statistical Me- 
chanics" Addison Wesley (Advanced Book Classics), 
Redwood City, 1991 

G.S.Setlur, M.T. Kuo and Y.C. Chang "Corrections 
to RPA Induced by Number Fluctuations", APS e- 
print, 1998 



G.S.Setlur an d Y.C. Chang, LANL preprint sond- 
mat/98103Qg| 



The pathological exception is for nonideal distribu- 
tions such as, (n k ) = for |k| < 0.5 kf and (n k ) = 1 
for 0.5 k f < |k| < 1.5 k f and (n k ) = for |k| > 
1.5 kf this is very nonideal but has zero fluctuation, 
but hopefully these are not physical. 



4 



